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Abstract 

We discuss the scaling behavior of the self-organized critical forest-fire model 
on large length scales. As indicated in earlier publications, the forest-fire 
model does not show conventional critical scaling, but has two qualitatively 
different types of fires that superimpose to give the effective exponents typ- 
ically measured in simulations. We show that this explains not only why 
the exponent characterizing the fire-size distribution changes with increas- 
ing correlation length, but allows also to predict its asymptotic value. We 
support our arguments by computer simulations of a coarse-grained model, 
by scaling arguments and by analyzing states that are created artificially by 
superimposing the two types of fires. 
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I. INTRODUCTION 



During the past years, systems which exhibit self-organized criticality (SOC) have at- 
tracted much attention, since they might explain part of the abundance of fractal structures 
in nature M. Their common features are slow driving or energy input and rare dissipation 
events which are instantaneous on the time scale of driving. In the stationary state, the 
size distribution of dissipation events obeys a power law, irrespective of initial conditions 
and without the need to fine-tune parameters. Examples for such systems are the sandpile 
model [Q], the self-organized critical forest-fire model (FFM) the earthquake model by 

Olami, Feder, and Christensen ||, and the Bak-Sneppen model @||. 

The study of SOC models is usually based on the assumption that the size distribution 
n(s) of dissipation events (avalanches, fires, earthquakes) shows the scaling behavior familiar 
from equilibrium critical systems, 

n(s) ~ s _T C(s/s max ) , (1) 

with a cutoff function C that is constant for small arguments and decays exponentially fast 
when the argument is considerably larger than 1. The cutoff cluster size s max is related to 
the correlation length £ via s max ~ £ D , with D being the fractal dimension of the dissipation 
events. (If the cutoff is set by the system size L, £ must be replaced with L.) This holds 
indeed for some self organized critical systems, like the Bak-Sneppen model, but it has been 
known for some time that it does not hold for the two-dimensional Abelian sandpile model 
||. Very recent work has shown that this violation of simple scaling in the sandpile model is 
due to the existence of multiple waves of topplings, and some features of the correct scaling 
behavior have been worked out |10|JTT[| . Violation of finite-size scaling is also seen in the 



above-mentioned earthquake model fT2|| . 



During recent years, evidence has accumulated that the two-dimensional SOC forest- 
fire model does not show simple scaling either. Instead, there are more than one diverging 



length scale [13], the behavior of the model for tree densities just above the critical density is 



2 



completely different from that of conventional critical systems [p^ , |15| , and finite-size scaling 
is violated |]IE|| . A scaling collapse based on Eq. (|1|) gives a good overlap of the tails of the 
distribution, but not so much of the first part, where the slope (i.e., the exponent r) seems 
to increase slightly with increasing correlation length (see, e.g., and Figure below). We 
have suggested [[I1J that all these features are due to the fact that there are two qualitatively 
different types of fires in the system: smaller, fractal fires that occur in regions of low tree 
density and burn a tree cluster that resembles a percolation cluster, and larger compact fires 
that burn a patch of a tree density above the percolation threshold. 

It is the purpose of this paper to show how these two types of fires add up to give the 
distributions typically seen in computer simulations, and to derive the asymptotic properties 
of the fire size distribution in the limit of very large correlation length. In particular, we 
will derive the asymptotic value of the exponent of the fire size distribution, toward which 
it should converge for sufficiently large correlation length. We support our arguments by 
three different types of evidence that all lead to the same conclusions: (i) Scaling arguments 
based on the superposition of the mentioned two types of fires; (ii) The fire size distributions 
resulting from the (artificial) superposition of the fires of patches of different tree density 
and different size; (iii) Computer simulations of a coarse-grained model that allow us to 
study numerically systems with much larger correlation length than has been possible so 
far. The outline of the paper is as follows: In the next section, we will derive scaling laws 
and analytical expressions for the fire size distribution resulting from the assumption that 
the forest-fire model is composed of patches of different size and different tree densities. 
Then, in section |T| we will show numerical data that result from the superposition of fires 
from artificially generated patches of different sizes and densities. In section [TV], we will 
present and study a coarse-grained forest-fire model where each lattice site stands for a 
group of several sites in the original model. Finally, we will discuss our findings. 

II. SCALING PROPERTIES OF THE SOC FOREST-FIRE MODEL 
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A. The definition of the model 



The SOC forest-fire model is usually studied on a square lattice with L 2 sites. Each 
site is either occupied by a tree, or it is empty. At each time step, all sites are updated 
in parallel according to the following rules: (i) An empty sites becomes occupied by a tree 
with probability p; (ii) A tree is struck by lightning with probability /. This tree and the 
whole cluster of trees connected to it (by nearest-neighbor coupling) burn down and become 
empty sites. 

B. Two types of fires determine the dynamics of the SOC forest fire model 

As long as p and / are so small that fires do not interfere with each other or with tree 
growth, the stationary behavior of the model depends only on the ratio f /p, but not on the 
two parameters separately. After some time, the system reaches a stationary state with a 
mean tree density p and a mean fire size s. A snapshot of the 2d system in the stationary 
state is shown in Figure [1]. One can see that it consists of patches of different tree density 
and different sizes. Some of the patches have a high tree density, and if they are struck by 
lightning, the entire patch burns down, with only few trees being left. After the fire, the 
tree density p pa tch(^) of the patch grows again according to 

PpatchO) = P(l - Ppatch(O) ( 2 ) 

until it is hit by the next lightning stroke. This mechanism of growth and burning down of 
forest clusters produces the patchy structure seen in Figure [j], which is characterized by the 
following properties. 

• The patches are almost homogenously covered with trees. This is because a fire that 
burns a patch usually leaves only a few trees behind. (We found that the local tree 
density within a patch immediately after a fire is typically 0.078 for the 2d square 
lattice and 0.062 for the 2d triangular lattice). Thus the random tree growth leads to 
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a uniform tree density within the patch. The size distribution of tree clusters within a 
patch is therefore similar to the size distribution of clusters in a percolation system of 
the same density. (For an introduction to percolation theory, see ||17|| .) The fact that 
several patches contain larger dense tree clusters indicates that some fires leave behind 
small clusters that can become seeds of new patches. This process of birth and death 
of patches, which happens on a slower time scale, does not affect our main argument. 

• It can be assumed that the distribution of patch sizes is independent of their tree 
density, because a fire in most cases hardly changes the size of a patch. Thus, the size 
of most patches is the same for high and low tree density. 

• Some patches have a tree density above the percolation threshold. These patches 
contain a spanning cluster that is compact, i.e., has the fractal dimension 2. When such 
a patch is struck by lightning, a compact fire occurs that burns the spanning cluster. 
If lightning strikes a patch of low tree density (below the percolation threshold), only 
a small, fractal cluster of trees burns down, and only part of the patch is affected by 
the fire. If the mean fire size is large (i.e. if f/p is small), most trees burn down during 
large fires, and most of the empty sites are created during the large compact fires, 
resulting in the above-mentioned low local tree density immediately after a fire. 

• The size of the largest patch diverges for f/p — > 0, suggesting that the system is close 
to a critical point and can be characterized by power laws. Several such power laws 
will be mentioned further below. 

The two mentioned types of fires add up to give the distributions typically seen in 
computer simulations, and explain the unconventional behavior of the forest fire model 
mentioned in the previous section. The left part of the fire size distribution of the forest 
fire model (see Figure 0b) is mainly due to fires burning fractal percolation clusters, and the 
cutoff part is due to large compact fires. In contrast, in conventional critical systems the 
power-law part and the cutoff part are due to the same type of critical fluctuations. Since 
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the two parts of the distribution become clearly separated only for very large £, as we will 
show below, the asymptotic exponent of the fire-size distribution is not visible in present-day 
computer simulations. 



C. The scaling behavior 

In the stationary state, the mean number of growing trees must equal the mean number 
of trees burning down, leading to 0] 

S-^A. (3, 

s diverges according to a power law in the limit f/p—>0, implying that the size s max of the 
largest fires also diverges, and with it the correlation length £, which we define to be the 
radius of the largest fires. The cutoff fire be expected to scale as 

Vax ~ (fM^, (4) 

with an exponent A 0-0, which has a value close to 1.1. This leads to £ ~ (f /p)~ X ^ 2 since 
the fractal dimension of the large fires is 2. (In contrast, earlier work was based on the 
assumption that large and small fires have the same fractal dimension, which was found to 
be D ~ 1.96 ||, with some authors not ruling out the value D = 2 @H.) 

Let s be the number of trees burnt during a fire, and n(s) the size distribution of fires 
in the system, normalized such that / n(s)ds = p. Since each tree is struck by lightning 
with the same probability, the size distribution of tree clusters is proportional to n(s) / s. We 
write 

n(s) = ni(s) + n 2 (s), (5) 

with rii(s) being the contribution from the smaller, fractal fires, and n 2 (s) being the con- 
tribution from the compact fires that burn an entire patch. If there was no qualitative 
difference between the two types of fires, the scaling law Eq. (JJ) would hold, from which one 
could derive the scaling relation 
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A = — T , (6) 

(This is obtained from the condition s = sn(s)ds/p and is given in many earlier publica- 
tions.) However, since we have two different types of fires, we cannot expect this relation to 
hold. Figure §o shows the size distribution of fires for different values of f/p. (The system 
size L was chosen large enough to avoid finite-size effects.) One can see that n(s) becomes 
steeper with decreasing f/p, and one can expect the slope to grow further until it reaches 
some limit value. The steepest slope occurring in this figure has the absolute value r « 1.3, 
which must be a lower bound to the asymptotic value of r. (The value r ~ 1.14 given in 
many earlier publications was obtained by taking some average value along n(s), which was 
less than its steepest slope.) A scaling collapse (see Figure |2]a) shows that the cutoff parts of 
the curves superimpose nicely, allowing to derive a value A ~ 1.1, which does not fit together 
with the scaling relation Eq. (|J). The same result for A was obtained by Pastor-Satorras and 
Vespignani |18[] using a moment analysis, confirming that the cutoff shows simple scaling 



behavior. 

Next, let us discuss the properties of rii(s). As f/p decreases, the first part of n(s) 
does not change any more. This indicates that rii(s) reaches an asymptotic form n*(s) as 
f/p—>0, with a cutoff that depends on f/p. We therefore write 

m(s) = ri*(s)C 1 (s/s maX)fracta i) = nl(s)Ci(s(f /p) Xl ), (7) 

introducing the cutoff function Ci(s/s m ax,fractai) for the distribution of the fractal fires and 
assuming that the maximum fractal fire size s max .f racta i scales with an exponent Ai. For 
sufficiently large s, n*(s) will reach an asymptotic power law with the "true" exponent 
t. We can estimate the value of r from the following argument: The large fractal fires 
stem from the percolation clusters in those patches that have a tree density close to the 
percolation threshold p pcTC . Thus the probability density of finding a cluster of size s is 
proportional to the probability that p pa t c h °f a large patch is large enough that percolation 
clusters of size s exist, multiplied by the probability density to find a cluster of size s in a 
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system at the percolation threshold. The fire size distribution is proportional to s times the 
cluster distribution, as we mentioned above. The probability to find a cluster of size s in a 
percolation system is determined by the size distribution of percolation clusters: 



n 



perc 



s) — s p Cp erc (s/s max p erc ) (8) 



with 

^max,perc (Pperc Ppatch) P (9) 

and <7perc = 91/36 ~ 2.528 and r perc = 187/91 ~ 2.055 In these patches close to 

the percolation threshold, the tree density increases with time approximately as p pa tch = 
p(l — Pp erc ). Therefore, the probability that the density is within a distance p perc — p pa tch of 
the percolation threshold is proportional to p perc — p pa tch, and the probability that a patch 
has tree clusters larger than s is proportional to s~ 1 ' <Jperc (see Eq. (§)). The probability that 
a fire of size s occurs is consequently proportional to 

— Tperc 1/ <-?perc ^ ^ 1.45 

This means 

t ~ 1.45. 

Next, let us estimate the cutoff exponent Ai. As we have seen above, the radius £ of 
the largest patches is proportional to (/ /p)~ A / 2 . The size of the largest fractal tree clusters 
is therefore proportional to (f /p)~ DpcicX ^ 2 , with the fractal dimension D peic of percolation 
clusters ~ 1.56, implying Ai 0.86. Beyond the cutoff size for fractal clusters, proportional 
to (f/p)~ Xl , the size distribution of fires must be dominated by the compact fires and 
therefore by the size distribution of patches. Since Ai < A, the "bump" (which is dominated 
by the compact fires) should span a larger fraction of the fire size distribution for smaller 
f/p. Figure |2|b shows that this is indeed the case. 

Assuming that the patch size distribution n pat ch('S) scales, too, we suggest a scaling form 

n patc h(s) ~ s b m ^s- b C 2 (s/s max ), (10) 



with s being the number of sites in a patch, and s max oc (/ /p)~ x being the area of the 
largest patch. Since most of the system is covered by large patches, b must be smaller 
than 2, requiring the factor in eq. fllQ|) in order to normalize / sra patc h(s)(is. The size 
distribution n 2 (s) for the compact clusters depends on ri pa tch(s), but the relation between 
the two is non trivial. The reason is the following: Assume that a patch is struck by 
lightning always when its density is so far above the percolation threshold that it burns 
down completely. In this case, the size distribution of the large fires would be proportional 
to s times the size distribution of the patches. However, in this case patches would never be 
destroyed. On the other hand, patches merge from time to time with neighbors, when the 
neighbor reaches a density above the percolation threshold before lightning strikes the patch 
with the higher density. In order to obtain a stationary patch size distribution, patches 
must therefore be destroyed from time to time. This can only happen if they are hit by 
lightning with a non vanishing probability as long as their density is sufficiently close to the 
percolation threshold, such that smaller dense clusters of trees are left behind by the fire 
that can develop into small new patches. For this reason the size distribution of the large 
fires is different from s times the size distribution of patches. This will be seen also in the 
next section. 



D. The exponent 5 

Additional support for the picture that the fire size distribution is the sum of two qual- 
itatively different contributions comes from the scaling behavior of the tree density. It has 
been known for a long time that the tree density approaches its critical value according to 

(Pc-P)~(f/P) 1/S , 

with 1/5 ~ 0.5 @-[5|]. (The most recent and probably most accurate value is 0.47 |18|.) If 



the fire size distribution obeyed the scaling law Eq. ([J), one would expect 5 to follow from 

Pc- P= / s~ T ds 



ma x 



Assuming that conventional scaling (see Eq. ([[])) holds and using Eq. ([|) and (||) this leads 
to 

With the apparent value of r around 1.14, this would result in a value of 1/5 much smaller 
than 0.47. With the asymptotic value r = 1.45 (see below), 1/5 would have to be much larger 
than 0.47. It has been pointed out recently that the observed value of 5 makes corrections to 
simple scaling necessary and a second contribution to n(s) has been suggested, which 
has a larger exponent r but the same cutoff as the main contribution, and which becomes 
negligible for sufficiently small f/p and sufficiently large s. In contrast to these authors, we 
argue that there occur not merely corrections to scaling, but that the scaling behavior of the 
SOC forest-fire model is fundamentally different from simple scaling. For this reason, there 
is no relation between the exponents 5 and r, since there is no single exponent r describing 
the entire fire size distribution. Whether there exists another relation between 5 and the 
fire size distributions ri\ and ri2, we do not know. 

III. SUPERPOSITION OF DIFFERENT TYPES OF FIRES 

In order to show that the fire size distribution seen in simulations can be indeed the 
result of the superposition of the two mentioned types of fires, and in order to confirm 
that the asymptotic value of r is 1.45, we superimposed the cluster size distributions of 2- 
dimensional lattices that were homogenously covered with trees, and that had tree densities 
between p a fter the fire = 0.078 and p max = 0.625, with weights derived from Eq. (H). This 



kind of superposition was sugested by S. Clar in [p]q| . The values for p a f ter the fire and p 



'max 



were measured for instance in ||15|1 . But in addition to the superposition of the different tree 



densities we also superimposed different lattice sizes I, distributed according to Eq. fllCf ) with 
/ = yfs and cutoff / max = £ = ^/s max . The lattices thus represent patches of different sizes 
and densities. In order to find the value of the exponent b, we performed superpositions for 
20 different values of b from 0.1 to 2. The results did not depend very much on b as long as b 
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was smaller than 1, suggesting that the value of b is in the interval (0,1), but not allowing us 
to fix it more precisely. The results are only reproducible when the statistics are sufficiently 
good. For this reason, we had to superimpose 10 4 or more systems. 

The results are shown in Figure 0a and Figure ^b. Please note that these figures give 
cluster size distributions n(s)/s and not fire size distributions n(s), i.e., the exponents are 
larger by 1. Figure |]a shows that the apparent exponent r « 1.14 typically found in 
simulations of the 2d forest fire model is reproduced by the superposition. The smaller slope 
for small s, and the bump followed by the cutoff, are reproduced as well. We have performed 
this superposition also for a triangular lattice. This lattice is most easily implemented 
by taking a square lattice and including next-nearest neighbor couplings along one of the 
diagonals in each unit cell. As for the square lattice, the range of tree densities was obtained 
from simulations of the SOC model, and was found to cover the interval [0.062,0.534]. One 
can see that the apparent exponent r is the same as for the square lattice, explaining the 
"universality" of this exponent with respect to a change of the lattice type found earlier || . 
(All other figures shown in this paper are for the square lattice only. ) Figure |3|b shows that 
the distributions of the two cluster types separate for larger correlation length £, and that 
the slope of the part of the curve that stems from fractal clusters tends to r = 1.45 as we 
calculated in section |J. A similar effect will be found in the coarse-grained model discussed 
in the next section (compare Figure ||a). 

Our results show also that the size distribution of the largest fires is related to the size 
distribution of patches in a nontrivial way, as mentioned in the previous section. If all large 
fires did burn complete patches, the bump of the fire size distribution would have a slope 
—b + 1, which is positive. Since this is not the case, many tree clusters must be contributing 
to n(s) that are large but do not cover the entire patch. This consideration should hold for 
any value of s max . We expect therefore the value of p max to decrease slightly with increasing 
s max , such that there is always a non vanishing contribution to n(s) of clusters that are large 
but do not cover the entire patch. 

We conclude that the superposition of homogeneous patches reproduces important fea- 
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tures of the SOC FFM. It is also an efficient way of studying the regime of large correlation 
length, which is not accessible to direct computer simulations. 

IV. A COARSE-GRAINED FOREST-FIRE MODEL 

A. Definition of the model 

In order to be able to study larger systems, we introduced a coarse-grained model where 
each site stands for a group of sites in the original model. The variable at each site of 
this coarse-grained model is the local tree density p S ite, ranging continuously from to 
1. The rules of our coarse-grained model are the following: (i) the density at all sites 
increases per time step by a small amount Ap site = p(l — p S i te ); (h) Lightning strikes each 
site with a probability /. If the density of this site is below the percolation threshold 
Pperc = 0.59, nothing happens. If the tree density on a site struck by lightning is above 
the percolation threshold p pe rc, this site and the entire cluster of sites above the percolation 
threshold connected to it burn down. The density on a site after a fire is a random number 
between and r. The parameter r takes short-range fluctuations in the density into account. 
The smaller r, the smaller the density fluctuations. Smaller r means consequently that the 
density at each site is the average of a larger number of sites in the original model. We 
therefore expect the coarse-grained model for small r to resemble the original FFM for large 

B. Properties of the coarse-grained forest fire model 

Although the coarse-grained model is not exactly the same as the original FFM, it shares 
many of its features. Figure |] shows a snapshot of the coarse-grained model for r = 0.1. 
The figure shows a patchy structure similar to the one in Figure |l| In many patches one 
can see sites of two different densities. This indicates that lightning often strikes a patch 
before all of its sites (which cover a range of densities of the width r) have a local density 
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above the percolation threshold, leaving behind some sites with a density just below the 
percolation threshold. If this happens several times within the same patch, one can expect 
the patch to be destroyed and replaced by a set of smaller patches. Such processes of birth 
and death of patches are not considered further in this paper, but of course they occur also 
in the original FFM, as can be seen in Figure [I]. As we have mentioned above, creation 
of new small patches must occur in order to balance merging and growth of patches in the 
stationary state. 

Next, let us consider the size distribution of fires in the coarse-grained model. Figure [|a 
shows our simulation results for r = 0.1 and different values of f/p. One can clearly see 
that the slope becomes steeper with decreasing f/p and appears to approach a limit slope. 
Figure |b shows the slopes d\ogn(s)/ds as function of s, indicating that the predicted limit 
value 1.45 is indeed correct. 

Figure |5]c shows a collapse of the cutoff parts of the curves, giving A ~ 1.1, just as in the 
original forest-fire model. We also performed a moment analysis of the fire-size distribution, 
giving the same result A ~ 1.1. 

Figure ^| shows the fire size distribution n(s) for three different values of r, and for the 
same f/p = 0.01. For smaller r, the slope becomes steeper and the cutoff bump becomes 
more pronounced, indicating that for smaller r the coarse-grained model resembles the orig- 
inal model on larger scales. For smaller r, the cutoff becomes larger. The reason is that 
for smaller r a site of the coarse-grained model corresponds to more sites of the original 
model. The same lightning probability / per site in the coarse-grained model corresponds 
to a smaller lightning probability in the original model when r is smaller. 

Let us now estimate how many sites z(r) of the original model correspond to a site in 
the coarse-grained model with parameter r. From Figure |2]a we find that 

W = A(f/p)- x 

with A ~ 30. Similarly, we have for the coarse-grained model 

s m , x = B(r)(f/p)- x . 
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From Figure [|c we estimate 5(0.1) ~ 100, and from the data shown in Figure |6| we then 
obtain 5(0.2) ~ 66 and 5(0.5) ~ 44. Now, s max sites in the coarse-grained model correspond 
to z(r)s max sites in the original model, and / in the coarse-grained model corresponds to a 
lightning probability f/z(r) per site in the original model. Therefore we have 

B(r)(f/ P y x = A(f/pz)- x /z, 

leading to 

z(r) = (B(r)/A) l/{x - l) ~ (B(r)/A) 10 , 

resulting in z(0.5) ~ 46, z(0.2) ~ 2650, and z(0.1) ~ 170000. The length scales of the 
coarse-grained model are reduced by factors of the order 7, 50, and 400 for these three r 
values, compared to the original model. 

A direct comparison of a fire size distribution of the original model and one of the 
coarse-grained model with r = 0.5 confirms these findings. We searched for two fire size 
distributions such that the ratios of their / values and the ratios of their s max values are 
similar. This ratio turned out to be around 45, as shown in Figure |7|, and in agreement 
with the finding of the previous paragraph. Figure [7], shows also that the shapes of the two 
fire-size distributions, while similar, are not identical. Identical shapes cannot be expected, 
since the coarse-grained model is not completely identical to the original model on larger 
scales. For instance, inhomogeneities arising in the original model within an area of size 
z(r), cannot occur in the coarse-grained model. This explains the difference in shape on 
small scales. On large scales, the difference in the shape of the cutoff is probably due to 
the fact that the process of slow destruction of large patches is slightly different in the two 
models. In both cases, lightning strokes hitting the patch when its density is only slightly 
above the percolation threshold make the patch more inhomogeneous. In the coarse-grained 
model, this leads to sites belonging to two widely different density intervals, as mentioned 
further above in the context of Figure |j. In the original model, this leads to a couple of 
smaller dense tree clusters being left behind by the fire, as can be concluded from Figure [IJ. 

To conclude this section, our coarse-grained model, while not being exactly equivalent 
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to the original model, shows the features expected for the original model on larger scales 
and confirms in particular the universality of the exponent A and our conjecture that the 
exponent r has an asymptotic value around 1.45. The detailed mechanism of birth and 
destruction of patches is somewhat different in the two models and leads to different shapes 
of the fire size distributions at small s and for the largest s. 



V. DISCUSSION 

In this paper, we have argued that the fire size distribution in the SOC FFM is the 
result of the superposition of two types of fires. The smaller ones are fractal percolation 
clusters, while the larger fires are compact and burn down a patch of a tree density above the 
percolation threshold. We supported this picture by a direct analysis of the model, by the 
artificial superposition of the two types of fires, and by the introduction of a coarse-grained 
FFM. 

One of our main results is that the asymptotic exponent for the fire size distribution 
is r ~ 1.45, and is visible only at length scales not accessible to present-day computer 
simulations. For values of the correlation length typically seen in computer simulations, the 
exponent r has an apparent value which is smaller, and which seems to be insensitive to 
the lattice type used in the simulations. Furthermore, we found that the cutoff exponent 
A has an universal value A ~ 1.1, which is measured in the original FFM as well as in the 
coarse-grained FFM for different simulation parameters. The robustness of this exponent is 
additionally supported by our earlier finding that the correlation length shows nice scaling 
behavior in a generalized model where trees can be immune to fire [JT9j . 



We could not find the precise form of the size distribution of patches, although we 
presented evidence that it should be characterized by an exponent b smaller than 1. The 
patch size distribution is the result of a slow and highly nontrivial process of birth and 
merging and destruction of patches. This process also determines the size distribution of 
the large fires, for which we could not give an analytical expression. 
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In contrast to the exponent r, we could not derive the value of the exponent A from 
analytical arguments. We cannot rule out that its true value is A = 1, and that there are 
logarithmic corrections which make it appear slightly larger than 1. 

As we have shown, several scaling relations familiar from conventional critical systems 
and in particular from percolation theory do not hold in the SOC FFM. Instead, the FFM is 
a new type of non equilibrium critical system that has no equivalent in equilibrium physics. 
It is characterized by different phenomena on different scales, and by a patchy structure 
indicating that neighboring sites tend to be synchronized by burning down during the same 
large fires. 

By introducing the coarse-grained model, we have shown that there exists an entire class 
of models that share the same main features of a patchy structure and two qualitatively 
different types of fires, the asymptotic exponent r ~ 1.45, and the cutoff exponent A ~ 1.1, 
while details like the precise shape of the cutoff and the precise mechanism of birth and 
destruction of patches may differ. 

Our results show that SOC in dissipative systems can be caused by mechanisms funda- 
mentally different from equilibrium critical phenomena. We expect that other dissipative 
SOC systems are driven to criticality by mechanisms similar to the ones found in the FFM. 
This applies in particular to the SOC earthquake model ||, where a patchy structure with 
partial synchronization of neighboring sites was also found pPf . Very recently, it was also 
found that this model contains two qualitatively different types of avalanches: those within 
a "patch", and those that enter it from outside and span the entire "patch" ||21|| . 
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FIGURES 
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FIG. 2. (a) Superposition of the fire size distributions n(s) of the SOC forest fire model for 



L = 1300 and f/p = 0.000118 (dashed line), 0.000169 (solid line), 0.000236 (solid line) and 
0.000394 (dotted line). In order to make the curves collapse, the vertical axis had to be scaled 
with {f /p)~ Xr v , using the effective exponent r av = 1.14, which can be interpreted as an average 
exponent over a certain range of s values, (b) Fire size distributions n(s) for L = 1300 and 
f/p = 0.0001183 (upper curve) and f/p = 0.001183 (lower curve). The upper curve was shifted 
vertically by a factor C = 10, in order to make the shapes of the two curves better visible. 
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FIG. 3. (a) Size distribution of tree clusters n(s)/sn(l) resulting from the superposition of 
lattices with l mSLX = 50 and b = 0.6 (solid line: square lattice; dotted line: triangular lattice), and a 
power law with the exponent r + 1 = 2.14 (dashed line), (b) The same for £ max = 2000 and b = 0.6, 
compared to a power law with the exponent r + 1 = 2.45. The tree densities we used cover the 
interval [0.078,0.625] for the square lattice and [0.062,0.534] for the triangular lattice. 
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L = 1000. The tree density of a site is represented by its grey shade, with larger densities being 
darker. 
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FIG. 5. (a) Fire size distribution n(s)/n(l) of the coarse-grained forest fire model for f/p = 



0.0031 (dotted curve), 0.005, 0.01, 0.0125, 0.025 and 0.05 (dashed curve), with the parameters 
r = 0.1 and L = 1000. (b) Slopes of the curves shown in (a), (c) Collapse of the fire size 
distributions of (a). In order to make the curves collapse, the vertical axis had to be scaled 
with (f/p)~ XrV , using the effective exponent r av = 1.25, which can be interpreted as an average 
exponent over a certain range of s values. 
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FIG. 6. Fire size distributions of the coarse grained model for r = 0.1 (dashed curve), r = 0.2 
(solid curve), and r = 0.5 (dotted curve). L = 1000 and //p=0.01 for all systems. 
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FIG. 7. Comparison of the fire distributions n(s) obtained for the original SOC forest fire model 
for / = 0.00000156, p = 0.01, and L = 800 (dashed line, C = 41, A = 1), and for the coarse grained 
model for r = 0.5, / = 0.00005, p = 0.01, and L = 1000 (solid line, C = 1, A = 44.9" 1 ). the ratio 
between the two cutoffs, A, is the same as the inverse ratio between the two lightning probabilities 
/■ 
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